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We consider the the van der Waals force involving excited atoms in general environments, consti¬ 
tuted by magnetodielectric bodies. We develop a dynamical approach studying the dynamics of the 
atoms and the field, mutually coupled. When only one atom is excited, our dynamical theory sug¬ 
gests that for large distances the van der Waals force acting on the ground-state atom is monotonic, 
while the force acting in the excited atom is spatially oscillating. We show how this latter force can 
be related to the known oscillating Casimir-Polder force on an excited atom near a (ground-state) 
body. Our force also reveals a population-induced dynamics: for times much larger that the atomic 
lifetime the atoms will decay to their ground-states leading to the van der Waals interaction between 
ground-state atoms. 

PACS numbers: 34.20.Cf, 31.30.jh, 42.50.Nn, 12.20.-m 


I. INTRODUCTION 

Casimir and van der Waals (vdW) forces are interac¬ 
tions between neutral macroscopic bodies or atoms aris¬ 
ing from the quantum fluctuations of both the electro¬ 
magnetic field and the atomic charges mm- They are re¬ 
sponsible for many characteristic phenomena in physics, 
chemistry and biology: the deviation from ideal-gas be¬ 
haviour in non-polar gases [3], latent heat of liquids, cap¬ 
illary attraction, physical absorption, and cell adhesion 
[3]. Dispersion interactions have even played an impor¬ 
tant role during the early stages of planet formation [5], 
and they are also supposed to have a fundamental role in 
selective long-distance biomolecular recognition [6]. Due 
to their strong distance-dependence, they become more 
and more important on the ever-decreasing scales of nan¬ 
otechnology, where they lead to the unwanted stiction 
of small mobile components [7j. In a series of ground¬ 
breaking experiments, the vdW force between an excited 
barium ion and a mirror have been measured to high 
precision IS El. Such experiments show an oscillating 
dependence of the vdW force on the ion. 

We will focus on the vdW force between two atoms, in 
excited states |n) and |Z). The interaction in this case is 
different from the interaction between two ground-state 
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atoms due to the possible exchange of a real photon be¬ 
tween the atoms. In the well-understood non-retarded 
regime, that is, for distances r much smaller than the 
wavelength of atomic electronic transitions, one finds 


Fnr (r) = - 


47r^eQr' 


E 

k,p 




E^-EA + EB - eP 


( 1 ) 


where, = r/r, are the matrix-elements of the 
dipole operator and E^ the energy relative to the state 
|fc). For downward transitions, Ep — Ep E^ — E^ can 
be negative and the resulting energy positive, yielding a 
repulsive interaction. Hence non-equilibrium situations 
can provide repulsive vdW interactions. 

The interaction at larger separations has been object 
of controversies. In a first group of works, it was pre¬ 
dicted that the magnitude of the retarded potential os¬ 
cillates as a function of interatomic distance [HHIll- In 
a later group of publications it was claimed that the re¬ 
tarded potential is non-oscillatory and proportional to 
1 /r^ [m ITMIT] . The conflicting results are due to subtle 
differences in treating divergent energy denominators in 
the photon propagators: the poles in the real axis can be 
avoided using the principal value prescription or adding 
infinitesimal factors in the energy denominators and this 
leads to different results. Both procedures are mathemat¬ 
ically correct, but they yield different physical results: a 
spatially oscillatory behaviour of the interaction in the 
first case and a monotonically decreasing behaviour in 
the second. 
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A group of recent works have used dynamical ap¬ 
proaches to address the problem. By an appropriate 
time-averaging procedure [18] or in the limit of vanishing 
atomic line widths [TU], a third result, for the vdW in¬ 
teraction on the excited atom, was found that oscillates 
in magnitude and sign. Note that an earlier approach 
based on time-dependent perturbation theory yields a 
non-oscillatory result for the force on the ground-state 
atom that is however valid only for times shorter than 
the lifetime of the excited state [20]. Similar consider¬ 
ations about timescales hold for the diagrammatic non¬ 
equilibrium description used in Ref. m- 

A very recent work claims that both results, the mono¬ 
tonic and the oscillating, are valid, but they describe dif¬ 
ferent physical processes |22|: the oscillating result is re¬ 
lated to a coherent exchange of excitation between the 
atoms, while the monotonic result is associated to a fast 
loss of excitation acquired from the initially excited atom. 
Another recent work finds that both forces can simultane¬ 
ously arise in a single set-up: the vdW interaction on the 
excited atom oscillates, in agreement with HHUHlIia, 
but the vdW force acting on the ground-state atom is 
monotonic |23j . This result would imply an apparent vi¬ 
olation of the action-reaction principle in excited systems 
in free space. However, it was shown that the momen¬ 
tum balance is restored when taking the photon emitted 
by the excited atom into account [24]. This emission be¬ 
ing asymmetric due to the presence of the ground-state 
atom, the emitted photon carries some average momen¬ 
tum, so that the difference between forces on the excited 
vs ground-state atoms can be interpreted as a photon 
recoil force. The situation is somewhat similar to the 
lateral Casimir-Polder force on an atom near a nanofi¬ 
bre, which is also associated with asymmetric emission 

m- 

In this paper, we study the van der Waals interac¬ 
tion involving excited atoms by means of a dynamical 
approach on the basis of the Markov approximation. We 
show that the damped internal atomic dynamics uniquely 
determines the oscillatory or monotonic behavior of the 
retarded interaction for excited atoms. In our dynamical 
model, the poles in the real axis are automatically shifted 
to the upper or lower part of the complex plane, and no 
ad hoc choice for the imaginary shifts in the denomina¬ 
tors is required. We will show that, when one atom is 
excited, the vdW force acting on the ground-state atom 
is monotonic and the vdW interaction of the excited atom 
is oscillating, in agreement with the most recent results in 
literature [2S] . Our dynamical approach is an al¬ 

ternative to the time-dependent perturbation theory |23j , 
where the behaviour of the force is determined via a time 
average over rapid oscillations on time scales of the or¬ 
der of atomic transition frequencies. Instead, our model 
allows us to study the decay-induced dynamics on larger 
time scales of the order of the excited-state lifetimes. It 
reveals that the force is governed by population-induced 
dynamics on these scales, where for times much larger 
than the lifetime of the intial atomic state the vdW force 


converges to that between ground-state atoms. In addi¬ 
tion, we are able to account for a general environment 
for the two atoms, via the classical Green tensor. 

The article is organised as follows. In Sect. [TTJ we 
present the basic formalism describing the coupled atom- 
field dynamics. It is used in Sect. |III| for calculating 
the force between two atoms in arbitrary excited initial 
states. In Sect. EYI we make the connection of our re¬ 
sult with the Casimir-Polder force between an excited 
atom and a body of arbitrary shape. Some conclusions 
are given in Sect. |Vj while in the Appendices, we present 
some of the more cumbersome details of our general ap¬ 
proach and our calculation. 


II. ATOM-FIELD DYNAMICS 

We consider the mutually coupled evolution of two 
atoms and the medium-assisted field. The field is pre¬ 
pared at zero temperature, and the atoms in generic in¬ 
ternal states. The dynamics of the atoms can be de¬ 
scribed with time-dependent flip operators, defined by 
Amra = where is an energy eigenstate, 

and similarly Bpg = 

In order to evaluate the force between the two atoms 
we must first solve the atom-field dynamics to obtain the 
flip operators and the field operators in the Heisenberg 
picture. The total Hamiltonian is the sum of three terms, 
the atomic and the field Hamiltonian and the interaction 
term in the multipolar coupling scheme within dipole ap¬ 
proximation: 

H =Ha + Hf + Haf 
Ha = E^Ann + Bnn 

n n 

oo 

Hf = J J dw/iwfl (r, uj) ■ h (r, w) 

X—e,m g 

HAF=-d^-E{rA)-d^-EirB) ( 2 ) 

where Ia (r, to) is the annihilation operator for the ele¬ 
mentary electric and magnetic excitations of the system 

m- 

Since the evolution of the whole system is unitary the 
commutator between two electric fields coincides with the 
commutator between free fields [55[ [H] 

[E (r, uj), E'^ (r', uj') 1 = ^^^ImG (r, r',uj) uj^S {uj — uj') 

(3) 

where G is the Green’s tensor of the electromagnetic field 
and E (r, uj) is the Fourier component of the electric field 
E(r) = /“dwE (r, w) -I- h.c., Heisenberg equations for 
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the coupled atom-field dynamics read 
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i9tE (r, w) = —iwE (r, w) 

ImG (r,r^, w) • + ImG (r,rB,a;) • d"® (4) 


where = 
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mn 5 



The electric field at the position of atom A consists of 
two terms: the radiation reaction and the field due to 
the other atom B. As shown in the literature El ED], 
the radiation reaction field gives rise to frequency shifts 
and spontaneous decay for atom A, see Figs. [^a) and 
(b). We thus renormalize the field by splitting off the 
radiation reaction 


= [id:;:^„ - (r;^ + r;^) /2] (i_) 

+ ^(K;^„-E 4 (r^,t)) (5) 

where m ^ n and the expectation value (...) is taken 
over atomic state and the field thermal state. E 4 {'^aA) 
is the sum of the free electric field and the source field of 
the atom B, the (second-order) Lamb-shifted atomic 
frequencies and F^ the decay rates. 


A 
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FIG. 1. Case a) is the zero order approximation: free field 
and free atom. Case b) is the next order approximation: the 
Lamb shift of an atom due to the emission and re-absorption 
of a photon. Case c) is the dispersion interaction between two 
atoms due to the exchange of two photons. 


III. VAN DER WAALS INTERACTION 
BETWEEN TWO EXCITED ATOMS 

We consider two atoms A and B that are initially pre¬ 
pared in excited energy eigenstates {ia), |*b) of the free 
atomic Hamiltonian (p(^(0) = 6niy^,pf{0) = Su^^). These 
initial states are not eigenstates of the total Hamiltonian 
and thus the atomic states evolve in time yielding a time- 
dependent vdW force (population-induced dynamics). 
As time progresses, the lower lying levels n < ^ is 

will become populated. 

To find the vdW force between on, say, atom A, we cal¬ 
culate the Lorentz force in electric-dipole approximation 
acting on A which is due to the field E 4 (r^,t) emitted 
by the other atom B: 

FA(rA,rB,t) = VA(d^ ■ E4 (r^,t) ), (6) 


where expectation value is taken over atomic and field 
states. 

For weak atom-field coupling, corresponding in an ex¬ 
pansion of the Hamiltonian in powers of the coupling 
strengths d, we can apply the Markov approximation to 
find (see Appendix): 


FAirA,rB,t) = 

n^iA B k,p 

00 00 

X J duj J dw'w^w'^Vyi {d;^fc • ImG (rA,rB,w) • d^; 

0 0 

16 

xdjp • ImG (rB,rA,a;') • d^„}^ — + h.c. (7) 

2 = 1 * 

where p^{t) = (A„„ (t) ), pf{t) = (Bu (t) ) represent the 
atomic populations of states \n) and |Z). 

The energy denominators Di are listed in Table 1. 
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TABLE 1. Energy denominators. In this table, repre¬ 
sents the transition frequency between the virtual state \k} 
and the excited state |n), while ojpi represents the transition 
frequency between the virtual state |p) and |1). Furthermore 

± i (r^ + r;^) / 2 , ± i (r^ + rf) /2 

and = ca ± ie with e infinitesimal factor. F is the atomic 
line-width. 


Due to our dynamical treatment of the atom-field cou¬ 
pling, the result explicitly depends on atomic damping 
constants or line widths, and also an infinitesimal damp¬ 
ing for the photon frequency ca. These factors uniquely 
ensure the convergence of time-integrals. 

For excited atoms the energy denominators can ex¬ 
hibit poles, for photon frequencies being resonant to the 
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atomic ones. According to time-independent perturba¬ 
tion theory these poles would be situated on the real- 
frequency axis with the mentioned resulting ambiguities. 
In our dynamical approach, with the inclusion of the 
atomic line widths, the poles are automatically shifted 
to the lower or upper part of the complex plane leading 
to unique resonant contributions. 

The total vdW force acting on A consists in two terms, 
a non-resonant contribution arising from virtual photons 
exchange, and a resonant contribution which corresponds 
to a possible emission of real photons by the excited 
atoms: 


Fa {rA, rB,t) = ^B,t) + (r^, r^, t) (8) 


In the limit of vanishing line-widths, the non-resonant 
contribution in an arbitrary magnetoelectric environment 
reads (see Appendix): 


F;?' {rA,rB,t) = ^ J d^5‘‘VATr{Q:^ (i^) 

0 

• G {ta, rs, iO -as (i^ ’ G (rs, r^, i^} (9) 

where we have defined the following polarizabilities of the 
initially excited atoms: 
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( 10 ) 


The resonant contribution reads: 


F"4 {rA,rB,t) = 

Mo E 

n^iA k<.n 

xdgfc-G(rA,rB,wgfe)-aB (wgfc) -G (fb, ta, • dg„} 

+ MoE^^fWE^^{K)' 

P<1 

X dg • G (rB,rA,a;g) • aA (wg) • G* (rA,rB,wg) • dg} 

( 11 ) 


atoms in free space: 
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( 12 ) 


where Xnk = and yi^ = rwg/c, e^. = r/r. When 

both atoms are excited, the monotonic and oscillating re¬ 
sults both contribute and can be attributed to different 
physical processes [22]: the oscillating result is related 
to a reversible exchange of excitation (“pendulation”) 
and the monotonic form with an effectively irreversible 
(Forster) excitation transfer. 

When only one atom is excited, the force acting on the 
excited atom is oscillating; on the other hand, the force 
acting on the ground-state atom is monotonic, coherently 
with the perturbative result in EDj. This implies a viola¬ 
tion of the action-reaction principle in excited systems in 
free space. The interaction is accompanied by the trans¬ 
fer of linear momentum to the electromagnetic vacuum; 
this momentum is ultimately released through directional 
spontaneous emission of the excited atom [24] . 

In Fig. Q, we show the vdW force acting on a ru¬ 
bidium atom and on a Cesium atom in free space, the 
Rubidium atom being in the excited state 5^Pi/2 and 
the Cesium atom in the ground-state 5^5'i/2 (see [21]); 
the force is represented for times much shorter than the 
atomic lifetime and much larger than the inverse of the 
atomic frequency, so that the populations of the states 
may be considered constant and the atomic dynamical 
self-dressing is not present. At large distances the reso¬ 
nant term dominates and the force on the excited atom 
shows Drexhage-type oscillations with an amplitude r“^. 
The force acting on the ground-state atom is monotonic. 
At small distances, we find a non-oscillating repulsive 
force for both atoms. 


F[N] 



For large distances the resonant contribution dominates 
over the non-resonant one. Two terms, one oscillating 
and one monotonic, are involved in the resonant contri¬ 
bution. Their behavior can be seen explicitly for isotropic 


FIG. 2. vdW interaction between one Cesium atom in 
the ground state (5^S'i/2) and an excited Rubidium atom 
(5 ^Pi/2 )- The thick line represents the force on Rubidium 
and the dashed one that on Cesium. 
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However our theory is more general because it includes 
the presence of general environments for the two atoms, 
like magnetodielectric bodies. Many differences arise in 
this more general case. Firstly the interaction can be 
described as a two-photon process, where the photons 
can be reflected by the body’s surface (see Fig. [^); this 
reflection is mathematically described in our formalism 
by the scattering Green tensor, which is known for many 
geometries and magnetodielectric properties. Secondly 
due to the presence of the additional body the action- 
reaction principle is also violated for ground-state atoms, 
with the interaction being accompanied by the transfer 
of linear momentum to the body. Lastly the total force 
acting on one molecule is not parallel to the interparticle 
separation vector. 



FIG. 3. Body-assisted vdW interaction: the exchanged pho¬ 
tons can be reflected by the body’s surface. 


than the atomic lifetime the force is repulsive and roughly 
one order of magnitude larger than the ground-state 
force. 

As stated above, the resonant force on the excited atom 
may be associated with photon recoil due to spontaneous 
emission. The fact that this force is stronger at small 
times can be understood from its ensemble-average na¬ 
ture: the probability of photon emission (and hence re¬ 
coil) is highest for small times where a large fraction of 
the ensemble atoms are still in their excited state. 


IV. COMPARISON TO CASIMIR-POLDER 
FORCE 

Let us compare our result with the experimentally 
observed single-atom Casimir-Polder force. If an ini¬ 
tially excited atom A is placed near a magnetodielectric 
body, the resonant contribution, associated with a possi¬ 
ble emission of a real photon, reads [28] : 

Fa (^a, t) = Ho (^) H 

n^iA k<.n 

X VARe{d(^fc • (rAAA.i^nk) ' ^kn] (13) 

where is the body’s scattering Green’s tensor. If the 
body is made up of ground-state atoms with polarizabil¬ 
ity O-B (w) and positions vb and number density 77(r), it 
can be expressed in terms of a leading-order Born expan¬ 
sion |5^ : 


We see that the resonant contribution vanishes for 
times much larger than the atomic lifetimes (F)^ oc 

Pn (0 = when the atoms have decayed to the 

ground-state. Fig. represents this population-induced 
dynamics for the force acting on the excited Rubidium at 
a given distance. We see that for time much larger than 

F(t)/F(0) 



FIG. 4. Population-induced dynamics for the vdW force act¬ 
ing on the excited Rubidium atom (thick line). The Rubid¬ 
ium atom is in the 5'^Pi/2 state, while the Cesium is in the 
ground-state. The distance between the atoms is r = lOnm. 

the atomic lifetime the force converges to the ground- 
state force, which is attractive. For times much smaller 


G^(rA,rA,w) d^rs ?7(rB)G° (r^, r^, w) 

• as (w) • G° (rB,rA,a;) -b ... (14) 

where G° is the free-space Green’s tensor. The substi¬ 
tution of this expansion into the single-atom Gasimir- 
Polder force leads to a resonant force 


Fa (^a, 0 = y 'q{vB)F\ (i^A, rs) 

= J dVB7?(rB)Mo X] 

n^lA k<. 7 l 

•G^ {rAAB,UJnk) - OLE (uj^k) (vBAA.iAl^k) -dfcn} 

(15) 


on the excited atom which is simply the sum over the 
(oscillating) resonant forces © on the excited atom due 
to the ground-state atoms constituting the body. Note 
that a monotonous force contribution is absent from the 
single-atom Gasimir-Polder force (15), as the atoms in 
the body are not excited. In the present combination of 
an excited atom interacting with a ground-state atom, 
we would expect the force on the body to contain a 
monotonous Gasimir-Polder force component. However, 
the force on the body is usually not considered in the 
context of Gasimir-Polder physics due to the strongly 
asymmetric mass ratio. 
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V. CONCLUSIONS AND OUTLOOK 

Our dynamical theory has allowed us to study the 
vdW force involving excited atoms in generic environ¬ 
ments. It is able to give a unique answer to the old 
puzzle whether the respective interaction is oscillating or 
monotonic, without recourse to ad hoc assumptions or 
prescriptions. 

When one atom is excited we have shown that the 
van der Waals force acting on the excited atom indeed 
shows Drexhage-type oscillations, while the force acting 
on the ground-state atom is monotonic. We have explic¬ 
itly demonstrated that the oscillating force is consistent 
with the respective Casimir-Polder force between an ex¬ 
cited atom and a ground-state body. On the contrary, the 
monotonic forces components cannot be deduced from 
the atom-body force in this way, because they act on the 
atoms inside the body whereas Casimir-Polder calcula¬ 
tions are usually restricted to calculating the force on the 
single atom in front of the body. 

The oscillating force on the excited atom could have 
profound implications on the spatial correlations of ex¬ 
cited atomic ensembles, in particular for Rydberg sys¬ 
tems. In addition, both the oscillating and monotonous 
force components are expected to arise in waveguides as 
recently studied in Refs. [33HSS1- K could be also interest¬ 
ing to generalize our model to include finite temperature, 
by chancing the fluctuation relations of the electromag¬ 
netic field, and to consider many-body vdW forces. 
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VI. APPENDIX 

A. Perturbative expansion of the force operator 

We consider the dynamics of an operator p{t), which 
is a superposition of operators 0„(t) with complex coef¬ 
ficients 


n 


and we introduce a time limit which acts only on the 
operators: 


(17) 

n 


The operators On{t) evolve dynamically according to the 
Heisenberg equations: 

dtOn {t) = ^ [dn (t) (18) 

where H is the total Hamiltonian. This equation can be 
integrated from the initial time to fo ^ given time t: 

t 

On (t) = On (to) + ^ y dil [On (h), H (tl) ] (19) 

^0 

This equation shows that the dynamical evolution of the 
operators p is: 


t 

p(t) = P(t)\t^to + y / (ii)] (20) 


The vdW force operator acting on the atom A, due to 
the presence of other atoms, is: 


FA(t) = V[dA • E4 (r,t)]^^^ 


( 21 ) 


where the field E^(rA,t) represents the total electric 
field, excluding the radiation reaction of atom A. Using 
Eq. (20) we find the dynamical equation for the force: 


F^(t) =F^(t) 


t—^to 






( 22 ) 

This equation can be reiterated considering now the dy¬ 


namics of the commutator 




which is 


a superposition of operators at the time H. Therefore, 
for weak coupling, we can construct a perturbative ex¬ 
pansion Fyi(t) in terms of operators at the initial time 
to- 


Fa (t) = Fa (t) 


t—^to 


I 

(ti)]^^^^^ 


t ti 


to to 

[FA(t)l^^^,H(h)l^^^^,H(h) 


. t2—>-to 


+ ... (23) 


In our model the electric field and the flip-operators of 
the two atoms are the dynamical variables of the system. 
Two different time-scales are observed for the dynamical 
variables; there is a fast free dynamics and a much slower 
dynamics due to the interaction between the atoms and 
the field. For example the free evolution of the flip opera¬ 
tors is on time scales of = 10“^®s while the dynamics 
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due to the interaction is on time scales of 
We define new dynamical variables according the formu¬ 
las: 

E' (r, w, t) = e“^‘E (r, w, t) 

(24) 

where: 

(25) 

The new dynamical variables change on the time-scale of 
the interaction and have the following commutator with 
the total Hamiltonian (See Eq. |^[^: 


where the expectation value is taken over the atomic-|- 
field free state |/). This approximate solution to the 
coupled dynamics is equivalent to an iterative use of the 
atom-field equations, and it is valid for weak coupling 
between atoms and field. 

The expectation value on free atomic and field states 
can be easily performed, since after the limit t 2 N-i —>■ 
the resulting operators are evaluated at the same initial 
time toi which represents the time at which the electric 
field and the atoms are uncoupled. 


TT 


X (t) w^ImG (r, rs, w) (t ), 


(0 : H (t)J = - J dujf^„ i-t) 
0 

X (t) ■ E'4 (r^, w, t) 

+e“^*E'^(rA,w,t) (t) 


(26) 


where: 


it) = i^'rnk (t) f^k it) ^tk 

k 

-Kn it) iL it) dL ) (27) 


K" 


and a normal ordering prescription is used. 

From Eqs. (26) we see that the commutator be¬ 


tween the Hamiltonian and a dynamical variable in¬ 
creases the number of electric dipole moments by one. 
Hence Eq. (23) represents a perturbative expansion of 


the force with the dipole element as perturbative pa¬ 
rameter. In particular the electric vdW iV-body force 
Fa (rA, ri...., tat-i) acting on A due to the other atoms, 
with positions ri....,r^r-i, contains 2N electric dipole 
matrix elements; this force results from the application 
of 2N — 1 commutators: 


Fa (rA,ri....,r7v-i,0 


2N-1 


ti t2N-2 


J dti J dt2.: J dt2N-i 

to to to 


X I 


[FAit)l^^^,Hih)]^^^^^,Hit2) 


J t2^t3 


., H {t2N-2) 


, H it2N-l) 


t2N-2^t2N-l 


t2N-l^to 


( 28 ) 


B. Van der Waals interaction between two atoms 


We consider now the vdW interaction between two 
atoms. We suppose that the atomic states are incoherent 
superpositions of energy eigenstates |n^) and \l^') and 
the state of the field is the ground state. 


In normal ordering, the force operator (see Eq. 21) can 


be expressed in terms of the new dynamical variables: 

OO 

Ea it) = ^J dw/^„ it) e-“^* 

m,n g 

X it) ■ E;^ (r, w, t) + h.c. (29) 


The two-body vdW interaction, which contains four elec¬ 
tric dipole moments, involves three commutators: 


g i i 1 i 2 

FA(rA,rB,t)= J dti J dt2 J dt^ 

to to to 


,H(t2) .H(h) 

-I t2^t3 


- to—^to 


(30) 


With the help of Eqs. (26), the commutators can be 


evaluated. For example the application of one and two 
commutators gives: 






















ih 


oo oo 

[F^ (t-ii)v{(e--'‘^K;^„ (h) ■±'^irA,u;\h) + 

m,n p Q 


J r=rA 


+ 


,,L,,L 0 0 

(r^, w', h) ■ (ti) )d;^„ • ±'^ (r, w, ti) 

OO 

I dwa;2e-“^(*-‘i (t)(ti) (ti) (ii)d;^„ • Im G (r, r^, w) • + h.c. 

m,n,r,s q 


iMo 

TT 


, ih. 


oo oo 

'[FA{t)Y^,H{h)]^^^^^,Hih)]=-^ Y fdLo fdu:'f^At-h)fr^At2) 

m,n,r,s q q 

^ ^ •ImG(rA,rB,w')(ta) 

- (ia) df^ • Im G (rs, r^, w') • K;^„ (ti) d;^„ • (r, w, ta) 

dmn • Im G (r, rs, w) • d^sKs ih) } 


X vjw'^e 




r—VA 


+ e*‘^'*iE'^(rA,w',t2) •K;^„(ii) d;^„ 

'-‘dt^ 2 y| jB (ia) • (rA,a;',ia) 


Mo 

ttM 


OO OO 

E 

m,n,r,s q q 


dw / dw'e-''^^*- 




= E'i (r^, w', ta) • KY it2) )B'rs ih) + fL A) frs ih - h) iU ih) 
X (ta) • E^ (rs, c^', ta) + e'^'^^E'^ (r^, co', ta) • Kf, (ta) )^ 

X d™„ • ImG (r,rB,a;) • dfJ + h.c. 

) r=rA 


(31) 


The commutators between K-mn and the Hamiltonian After some algebra we obtain: 
have not been considered since they lead to higher or¬ 
der corrections in the electric dipole d"^ and d^. 

We then evaluate the last commutator and take the 
expectation value on the atomic and field states. The 
thermal expectation value over the free field variables can 
be performed with the help of the following fluctuation 
relations for zero-temperature mM- 

(r,a;,t) E(°)I (r',a;',t)^ 

= ^^ImG (r,r',a;)a;^(5(a; — w') (32) 


TT 
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F {rA,rB,t) = J J dw'wV^V^d;^;, • ImG (rA,rB,w)-d^ 

n,l k,p Q Q 

t ti t2 

xdg-ImG {rB,rA,u}') ■ d^n J j dt2 J dts 

to to to 

X (/„! {t - h) - ft {t - t2)) 

+ {t - h) - {t - h)) (/,f (^i - t^) - f^i {h - h)) 

+ {ft {t - ti) - ft {t - ti)) ih - h) - e-(‘-*3) jB 

+ e--(‘-*3) (e--(‘-*3) - e- {ft {t - t,) - ft {t - h)) f§ {h - ta)} + c.c. (33) 


where is now applied to both Green’s tensors (after 
exploiting their symmetry and introdu cing a factor 1 / 2 ). 
The function / was defined in Eq. (251 and p)/(t) = 
{Ann{t)) and pf{t) = represent the atomic 

populations of the states \n) and |?). We have consid¬ 
ered time-reversal symmetric systems where d^n is real 
(dmn = d„m), and reciprocal media (G"*" (r^i, r^, w) = 
G (rB,r^,w)). 


With the exception of resonant cavity-QED scenarios, 
we can assume the quantity w^ImG (r^, r. 4 , w) to be suf¬ 
ficiently flat and to not exhibit any narrow peaks in vicin¬ 
ity of any atomic frequency (weak coupling). For weak 
coupling, we may evaluate the time-integral by means of 
the Markov approximation, extending the lower limit of 
the time integral to to = —00. The resulting integrals are 
not converging. In order to force the convergence we add 
an infinitesimal factor to the frequency w, uj —>■ a; — ie, 
where e > 0. Note that the opposite sign convention 
for this infinitesimal factor would lead to divergent inte¬ 
grals. Time-integration leads to the energy denominators 
in Table 1 in the main text. 


The frequency denominators can be combined: 


I/D2 + I/D7 + I/D 


10 

= l/(a;(") - 

I/D3 + 1/De + 1/^11 

pi 


— l/{t ^ + uj'){u}'+ 0 


1 /Di + 1 /De = 


1 


(wl '>+uj'){ujt ') 


^ ^ + ^pl' 


1/D4 -I- 1/Ili2 — 




I i(-) / ^ / 4- / 

’ - ^kn ^ + ^pl 


1/^5 = 


1 


{uj( )-uj'){ujt ^) 

1 1 


,I , A(~) , ,(—) I , 

^ +^kn 


1/De = 


1 


(a;(-)+020(0^^1+^ 

1 1 


^ + ^kn ' ~ 


^kn 


( 34 ) 


which implies: 


2 = 1 ^ 
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where we have defined the following functions: 
/i(?)=- ^ 


,^(+) , .B(.+) 


+ W, 


pi 




^kn J 


1 


(‘ 




'^pl 


')(« 


B(-) 

Si 


/2 (0 = 




1 


^(+) , , B(+) 


+ W, 


pi 


(« 


— Lu' I 


kn 


1 Wi.„' ' + a;„,' ' I (s 


M-) 

'^kn 


B{-: 

%i 


) 

A(-Y^ 


^kn 


I 

1 ^ 

\ 1 

- (36) 





= 

^pi ± 


1 ^(i) 

and a;fe„ 

i (r^ + rf) /2, a;(±) = w ± ie. 


For the first term in Eq. (35) we integrate over w and 
for the second term we integrate over w'. We use the 
identity ImG = (G — G*) /2i and the Schwarz reflection 
principle for the Green tensor: 


OO 

/ dw'w'^ ( ——^ ^ ) ImG (rA, rs, w') 

J yw'+w'- > uj'— oj^ > J 


= 2i / 


uj' + UJ^~'> 


+ 


io' — Ul^ 


G (rA,rB,a;') 

(37) 


The Green’s tensor is analytic in the upper half of the 
complex plane, including the real axis, and it is also finite 
at the origin. We close the path with an infinitely large 
half-circle in the upper complex half-plane and take the 
residuum inside the path. The integral along the infinite 
semi-circle vanishes for because: 

lim G(rA,rB,a;)|^ ^ =0 (38) 

We thus find: 

OO 

0 

~ 7ra;^G(rA,rB,-a;) (39) 

The total force can be expressed as sum of two terms: 
Fa (rA, rs, t) = F^ (rA, rs, t) + F^ (rA, rs, t) 


Fa irA,rB,t) =^J2Pnit)pfit)J2 J 

n^l k,p g 

X VAlm|(d(^fc • G(rA,rB,a;) • gi (a;)| 

OO 

Fa irA,rB,t) =^^p^(t)pf{t)Y^ j dww^ 

n,Z k,p Q 

X VA||d;^fc •G(rA,rB,w) •d^|% 2 (w)| (40) 

where: 

5i (w) = fi (w) -b /2 (w) = 

1 / 1 


U! +U}. 


B{+) 

'pi 


U1 + U1 


A{+) 

kn 


, M-) M-) 


,^(+) I , B{+) 

'^kn + ^pl 


A(-) .,B{-) 

'\n + ^pl 


A{+) 

^kn 


+ 


to + to. 


B{+) 

pi 


4n ^ “ ^kn , 


(41) 


and: 


g 2 (w) = Im [/i (w) -b /a (w)] 
= 2Re ^ ^ 


24( + ) . 24( + ) 

L^^kn ^ku 


Im- 


to + to. 


B{-) 

'pi 


(42) 


We consider then the limiting case of vanishing line- 
widths: 


eA = {y + Tt) /2^0+ 

eB=(rf+ r«)/2^o+ 


-pj,~ ■ - (43) 

In this limit the function gi can be simplified: 

^ 4 (w - a;^„) (w-b w^„) 

hm gi (w) = 7 ^- 

((w + y) + e\j (w - y + iea) 


fw -b -b topi 


+^pi +) ( y.+y 


0(‘ 


'kn ^ '^pl 


(44) 


Using the property = T’y =p iTrd (x), where V is the 
principal value, we can also simplify g 2 : 


g 2 {to) = 2 -k Re 


A{+) , A(-S) 

L^kn ‘^kn 


d (w - lofp) 


(45) 

With these results, after performing a Wick rotation 
on the imaginary axis we find the following non resonant 








































11 


and resonant contributions to the F^: 

OO 

F\ {rA,rB,t) = ^J2Pnit)pF{t)J2 J 

n,l k,p g 

X G(r.,re.it) dgr} 

+ aSpxWpPwlZ 

n,l k,p 

X V^lResi [gi {uj)uj^ (d;^^ • G {rA,rB,uj) ■ d®;)^] 

- Res 2 [gt i-uj) (d;^^ • G (r^, r^, w) • d^;)^ ] | (46) 

where Resi indicates the sum of the residues in the first 


quadrant, and Res 2 the sum of the residues in the second 
quadrant. Similarly, reduces to 

F\ {VA, rB,t)= E 

l P<1 

x{dip • G {rB,rA,ujPp) ■ ola {uJ^p) ■ G* {ya,yb,^^Fp) ' dpj 

(47) 

where a a is the polarizability of the excited atom A. The 
sum of FJj and F\ gives the non-resonant and resonant 
contributions of the total vdW force, see Eqs. and 

m. 
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